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Neutrino Factories and Muon Colliders' cooling lattices require both high gradient rf and strong 
focusing solenoids. Experiments have shown that there may be serious problems operating rf in 
the required magnetic fields. The use of high pressure gas to avoid these problems is discussed, 
including possible loss problems from electron and ion production by the passage of an ionizing 
beam. It is also noted that high pressure gas cannot be used in later stages of cooling for a muon 
collider. Experimental observations using vacuum rf cavities in magnetic fields are discussed, current 
pubUshed models of breakdown with and without magnetic fields are summarized, and some of their 
predictions compared with observations. 

A new theory of magnetic field dependent breakdown is presented. It is proposed that electrons 
emitted by field emission on asperities on one side of a cavity are focused by the magnetic field 
to the other side where they melt the cavity surface in small spots. Metal is then electrostatically 
drawn from the molten spots, becomes vaporized and ionized by field emission from the remaining 
damage and cause breakdown. The theory is fitted to existing 805 MHz data and predictions are 
made for performance at 201 MHz. The model predicts breakdown gradients significantly below 
those specified for either the International Scoping Study (ISS) [1] Neutrino Factory or a Muon 
Collider [2]. 

Possible solutions to these problems are discussed, including designs for 'magnetically insulated 
rf ' in which the cavity walls are designed to be parallel to a chosen magnetic field contour line and 
consequently damage from field emission is suppressed. An experimental program to study these 
problems and their possible solution is outlined. 



PACS numbers: 29.20.-c, 29.25.-7, 29.27.-a,52.59.-f,79.70.+q 



I. INTRODUCTION 

A. The use of rf in Magnetic Fields for 
Neutrino Factories and Muon Colliders 

Low frequency (330-200 MHz) rf is needed for 
phase rotation and early cooling in the currently pro- 
posed Neutrino Factory [1] and Muon Collider [2] de- 
signs. The magnitude of the required magnetic fields 
are of the order of 1.75 T for the phase rotation and 
around 3 T in the early cooling lattices. The rf gra- 
dients specified are between 12 and 15 MV/m [3]. 

Experimental data exists on the operation of 
805 MHz vacuum rf in magnetic fields. These were 
obtained with two very different cavity types: 1) a 
multi-cell cavity with open irises [4], and 2) a single 
'pillbox' cavity with irises closed with Cu plates or 
Be windows [5], [6]. Both showed significant prob- 
lems and will be discussed in this article. 

There is very little data on operations of 201 MHz 
vacuum cavities in significant magnetic fields, al- 
though experiments are underway with local fields 
of about 1 T on a limited part of the cavity. Tests 
are planned for the 201 MHz cavity in fields of the 
order of 3 T and geometries that can be made close 
to those in the machine designs. [32] 

There is also data on the DC operation of a test 
cavity in high pressure hydrogen gas that showed no 



magnetic field dependence [7]. However, there may 
be other problems arising when an ionizing beam 
passes through such a gas filled cavity; in addition, 
such gas filled cavities cannot be used in the later 
stages of cooling for a Muon Collider because the 
Coulomb scattering in the gas would cause too much 
emittance growth. 

We first discuss the gas filled cavity experiments 
and the possible problems with their use with ion- 
izing beams; then we consider the experiments and 
models of vacuum cavities ( no gas in them) in mag- 
netic fields. 



B. 805 MHz High Pressure Gas filled Test 
Cavity Experiments 

A simple test cavity (Fig. 1) has been operated 
with hydrogen at different pressures and with 'but- 
tons' made of different materials that define the 
small high gradient gap [7]. It was found that at 
lower pressures, the breakdown gradient follows the 
Paschen prediction, but at higher pressures the gra- 
dient is limited to values (50 MV/m for Cu) close to 
those observed (53 ^ in the open multi-cell cav- 
ity discussed in Sec. ID) in vacuum cavities after 
conditioning. This similarity in the observed max- 
imum gradients suggests that the initiating mecha- 
nism for breakdown in the high pressure gas, and 
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FIG. 1: (Color) Schematic of a test cavity. 



vacuum cavities is the same. Models for this phe- 
nomena will be discussed in Sec. II. No change in 
breakdown was observed in the presence of an ex- 
ternal magnetic field of 3 T. This is as would be ex- 
pected from the mechanism for vacuum breakdown 
with magnetic fields that is described in Sec. IE. 



C. Theoretical Expectation for the Effects of 
an Ionizing Beam passing through a High 
Pressure Gas filled Cavity 

Tollestrup [8] has studied the likely effects of a 
muon beam passing through gas filled rf cavities. It 
is concluded that the lifetime of the electrons and 
ions produced by the ionization of those beams are 
long compared with the likely duration of the muon 
beams. It is further concluded that such electrons 
in the cavity will be driven backwards and forwards 
as the rf voltage oscillates, and that this will lead 
to heating of the gas and loss of the rf energy. The 
final Q of three different cavities exposed to 10^^ 
muons are estimated and given in Table I. The final 

TABLE I: Final Q of three cavities. 



Frequency (MHz) 


400 


800 


1600 


rf gradient (MV/m) 


16 


16 


16 


Final Q 


647 


325 


163 



Collider Collaboration (NFMCC) 2008 meet- 
ing [9], [10]. 2 

• The phase rotation and initial cooling for the 
high emittance muon collider [10]. 

Table II gives the needed parameters of the three 
cases. For the two collider cases, the numbers of 
initial muons are estimated from the final number 
of muons, multiplied by two for the 2 charges, and 
divided by the quoted transmission. This is a some- 
what optimistic assumption since the phase rotation 
and initial cooling will also see unwanted kaons, very 
high energy pions, protons and knock-ons from neu- 
trons, that will add to the losses. For the neutrino 
factory case the total number of muons is derived 
from the number of protons per pulse at 10 GeV 
together with MARS [11] simulated total pion pro- 
duction at this energy. 

We assume that the final Q of the cavity is given 

by 



(1) 



where / is the cavity frequency, N is the number of 
muons per pulse and the subscripted values are from 
the 400 MHz case in Ref. 8. 

Since Q tells us the energy loss per period of the 
rf, the Q = 22 for the low emittance case means 
that little rf is left after only three cycles. Yet for 
the phase rotation schemes under study, the muons 
are spread over 10 or more cycles prior to bunch- 
ing and rotation. Further down the cooling channel 
the numbers of muons are less, but until and unless 
the bunches are merged, the rf must remain effective 
over 10 or more cycles, and the required Q must be 
well above 60 [33]. In any case, it must be noted 
that high pressure gas cannot be used in the later 
cooling stages where the emittance is very small. In 
these cases the Courant Snyder (3± must be very 
small where any material is introduced. This can 
be achieved in local areas, or inside very high field 
solenoids, but not over the lengthy rf systems. Fi- 
nal cooling for a muon collider will thus inevitably 
require vacuum acceleration near strong fields. 



D. Experimental Operations of a Multi-cell 
805 MHz Open Cavity in approximately Axial 
Magnetic Fields 



Q for other cases can be approximately estimated by 
scaling from these numbers. This is done for three 
cases: 

• The phase rotation and initial cooling for the 
ISS Neutrino Factory [1] 

• The phase rotation and initial cooling for the 
low emittance muon collider parameters as 
presented at The Neutrino Factory & Muon 



An ionization cooling lattice had been designed 
that employed liquid hydrogen absorbers inside high 
field solenoids separated by lower field matching 
solenoids surrounding 6 cell, open iris, rf cavities. 
Such lattices used early in a cooling sequence, would 
use 201 MHz, but higher frequency and high mag- 
netic fields would be used for later cooling. 

To test the system a 6 cell cavity (Fig. 2) was de- 
signed with iris apertures tailored to fit the beam 



TABLE II: Parameters for a Neutrino Factory and low and high emittance Muon Colhder. 



Frequency (MHz) 

rf gradient (MV/m) 

Repetition rate (Hz) 

Wall power (GW) 

10 GeV protons on target (10^^) 

Final muons per charge (10^^) 

Cooling transmission (%) 

Muons per proton 

Muons before cooling, N (10^^) 

Final Q 



u Factory 

200 
16 
50 
4 
51 



0.33 
17 
7.6 



Muon ColUder 
Low emittance High emittance 



200 
16 
65 
3.6 
(35) 
10 X 0.1 
30 
(0.17) 
6 

22 



200 
16 
12 
3.2 
(170) 
2 
7 

(0.33) 
57 
2.3 



^ loH 






-80 



-40 



Magnet Center 

\ 



-20 



-10 

z (cm) 



20 



Window 
Damage 



40 



FIG. 2: (Color) Schematic of a multi-cell cavity first tested in the lab G magnet at Fermilab. The lines shown are 
the magnetic field lines extending from two irises to the end window. 



profile in the matching fields [4]. The cavity was 
then tested with and without magnetic fields. The 
available fields were not those in the matching de- 
sign, and the cavity could not be positioned symmet- 
rically in the magnet, so the fields used in the test 
were asymmetrical. The maximum achieved surface 
gradients were quite high (^ 53 MV/m), compara- 
ble with those in the gas filled cavity, and were not 
strongly dependent on the magnetic fields (Fig. 3a). 
However, both X-rays and dark currents were greatly 
increased with the magnetic field on. After some 
time, vacuum was lost, and it was found that the 
end vacuum window was so damaged as to cause 
the vacuum leak. 

It was found that the location of the window dam- 
age corresponded to a focused dark current coming 
from one of the high field irises. In addition radiation 
damage patterns were observed showing that many 
beamlets were focused by the magnetic fields onto 
the end window. There was no indication that dark 
currents from one high gradient iris were focused on 
to another, which may explain why the maximum 
achieved gradients remained high (see Sec. IIIB). 



E. Pillbox Cavity Breakdown in Magnetic 
Fields 

In a linac with open irises, the peak surface fields 
are typically a factor of two higher than the average 



accelerating gradient. Subsequent to the design of 
the multi-cell open iris cavity, it was realized that 
with muons one could introduce thin Be windows 
at each iris and obtain accelerating gradients much 
closer to the maximum surface fields. In addition, 
the resulting 'pillbox' cavities give more acceleration 
for a given rf power. A test 'pillbox' cavity was thus 
designed (Fig. 3b) and tested in a magnetic field. 
This time the cavity was mounted in the center of the 
magnet and could thus be tested with symmetrical 
fields. It was also tested with one of the magnets 
coils unpowered to give asymmetrical fields and with 
the coils powered in opposite directions to study the 
effects with these different field geometries. 



Without field the maximum achieved gradients [6] 
were somewhat lower than for the multi-cell cavity, 
but this was, at least to some extent, the result of 
more conservative operation after the severe dam- 
age seen in the earlier cavity. With magnetic fields, 
the maximum gradients were found to be strongly 
dependent on field (Fig. 3a). It was also found that 
over time its performance deteriorated. Examina- 
tion of the inside of the cavity showed severe pitting 
on the irises. The Be windows themselves showed 
no visible damage, but there was a spray of Cu over 
their surface and Cu powder in the bottom [12]. 
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FIG. 3: (Color) Left: a) breakdown gradients vs axial magnetic field. Right: b) schematic of a pillbox cavity. (Data 
from Ref. 5) 



II. BREAKDOWN MODELS WITHOUT 
MAGNETIC FIELDS 

It is assumed in all models that breakdown is ini- 
tiated at 'asperities', where the local electric fields 
is higher, by a factor, ^pN introduced by Fowler- 
Nordheim [13]. The average values of these factors 
can be determined by observing the electron cur- 
rents (dark current) emitted by the sum of many 
asperities, each of which has a specific value of /3fn- 
The field emitted average electron current density 
Jf{^) for a surface field and local field 

£^iocai = PfnE is given by [14] 



Ai . lb 



C0 



1.5 



local 



where (j) is the material work function in (eV) 
4.5 eV for Cu) and 



C = 6.53 X lO^(eV) 



-1.5 



(2) 



(3) 



In vacuum cavities with a thin window, one can 
measure some fraction of all field emitted electrons 
and observe its field dependence. From this depen- 
dence, given an assumed work function 0, one can 
then extract an average value of (^^fn) 

It is reasonable to assume that breakdown oc- 
curs where the local field, and thus (3ip^E is max- 
imum, and thus higher by a factor a than the aver- 
age value determined from the gradient dependence 
of the dark current from many asperities, so 



local 



(/Sfn) {E) 



(4) 



where a>l depends on the probability distribution 
of ^FN. 



Observed breakdown gradients are found to de- 
pend on frequency [15] (oc vT)? pulse length, and 
cavity dimensions, but it has been found [16] that, 
over a range of frequencies from DC to a few GHz, 
and for differing pulse lengths, cavity dimensions 
and in waveguides, the values of ^^^^ fall in a rel- 
atively narrow range around 7 GV/m. Breakdown 
thus appears to be related to the local electric fields 
at asperities (points, cracks or other causes of the 
field enhancement), or to the field emitted currents 
that are strongly dependent on these fields. Several 
plausible but speculative mechanisms for the initia- 
tion of rf breakdown have been proposed, which will 
be discussed next. 



A. Breakdown Models 

1. Mechanical fracture model 

The mechanism, as discussed in Ref. 17, assumes 
the following sequence of events: 

1) The surface contains asperities at the top of 
which the local field is given by the average field 
multiplied by a Fowler-Nordheim field enhancement 
factor y^FN- 2) The outward electrostatic tension 
Fs = y(/^fn£^)^ is equal in magnitude to the en- 
ergy density of the field, hence it is proportional to 
the square of the field (7 GV/m would induce a ten- 
sion of 300 MPa), breaks off the tip and the small 
piece now moves away from the remainder of the as- 
perity. 3) The piece is bombarded by field emitted 
electrons from the remaining asperity and becomes 
vaporized and ionized. 4) Following this formation 
of a local plasma other mechanisms cause the plasma 
to spread, or other mechanisms short out the cavity 
leading to breakdown. 



Within this model the breakdown occurs when the 
electrostatic outward tension at the asperity equals 
the tensile strength of the material. For the local 
field ^^locai = <^(/^fn) {E), the average field at the 
surface of the asperity is given by 



{E}<x 



(/?fn) a 



(5) 



where T is the tensile strength of the material. 
There is relatively little data on breakdown, includ- 
ing y^FN, of materials of significantly different ten- 
sile strength. Figure 4a shows the observed depen- 
dencies in an experiment [18] at 11 GHz in a spe- 
cial tapered rf waveguide. The error bars on ten- 
sile strengths indicate the range of quoted values, 
depending on material treatment, with the highest 
strengths given for fine drawn wire. It is reasonable 
to assume that on the nanometer scale of an asperity, 
the material strength is of the order of that the high- 
est strengths observed. The plotted curve represents 
the theoretic predictions assuming a = 1. However, 
independently of the value of a, the theoretical ex- 
pectation is not observed in the experimental de- 
pendency; indeed, the higher strength stainless steel 
has a lower minimum local breakdown field than ex- 
pected. 



2. Ohmic heating model 

5 

It has been suggested [14] that, when the field 
emission current density is sufficiently high, break- 
down is initiated by ohmic heating that melts the 
tip of an asperity. Once liquefied electrostatic forces 
would pull the molten material away, just as the bro- 
ken piece of the asperity was pulled away in the first 
model. This molten material, as it lifts from the re- 
mains of the asperity, will be exposed to field emis- 
sion from the remaining asperity left behind and it 
will be further heated, vaporized and ionized to form 
a plasma. For submicron asperities, the time con- 
stant for achieving a steady thermal state is only 
of the order of a nanosecond, so the temperatures 
reached depend only on the geometry, electrical re- 
sistivity, thermal conductivity, and current densities 
at the tip. Assuming the asperity to have a conical 
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FIG. 4: (Color) Top: a) minimum local rf gradient at 
an asperity over a vs tensile strength of materials of an 
11 GHz waveguide; the error bars on tensile strengths in- 
dicate the range of quoted values, depending on material 
treatment. Bottom: b) local rf gradients vs the product 
of the melting temperature Tm times the thermal con- 
ductivity K divided by the electrical resistivity p for an 
11 GHz vacuum waveguide. (Data from Refs. 7, 18) 
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5: (Color) Schematic for asperity heating calcula- 



shape (see Fig. 5), with solid angle and emitting 
area A = h^Q^ then, given the electrical resistivity 
p and current density jo , the heat Q fiowing back to 
the base of the asperity is 



Q{h) = [ 



ho h 



(6) 



The temperature difference between the tip and the 
base AT, as a function of the thermal conductivity 
i^, assuming hi ho is 



AT = 



phi 
Jho 



Kh'^n 



dh 



I 



2hlKQ? J 
= ('^] (7) 



The field emission current density is approximately 
proportional to the local field £^iocai to the tenth 
power [4], so for a fixed emission area A and cone 
angle 1], the field needed to melt the material is pro- 

portional to £^iocai oc ( — ^ j . This expression 

contains a number of approximations: the temper- 
ature dependence of the parameters is ignored, the 
asperity shape is assumed to be a cone of fixed angle, 
and the field emission is approximated by a power 



law, but the result should be qualitatively correct, 
and is plotted in Fig. 4b. It is seen that the ma- 
terial dependencies agree with the prediction within 
the errors, and, in particular, the lower gradients 
achieved with stainless steel are as predicted. 

3. Thermal runaway model 
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FIG. 6: (Color) Rates of heating from electrical resistiv- 
ity p, and cooling from thermal conductivity vs the 
effective temperature at which p and K are determined. 

The above calculation of ohmic heating was car- 
ried out under the assumption that the resistivity 
and thermal conductivity are independent of tem- 
perature, which they are not. Assuming that these 
expressions are approximately valid using resistivi- 
ties and thermal conductivities at an intermediate 
effective temperature Teffective, then one can look at 
the relative rates of heating and cooling as a func- 
tion of that effective temperature. Figure 6 shows 
the heating and cooling vs the effective temperature 
for pure Cu [19]. For a low current density jo giving 
a temperature less than a critical value (Teffective < 
^critical), then the cooling which is ex AT x K rises 
more rapidly than the heating (oc p{T))^ and a sta- 
ble temperature is possible. On the other hand at 
higher current densities, which leads to temperatures 
above Tcriticai, the rate of heating vs temperature 
rises faster than the cooling and the temperature 
will 'runaway'. For the data used, the critical effec- 
tive temperature is about 300°. The actual critical 
temperature at the tip will be somewhat higher. It 
is worth noting that repeated heating to tempera- 
tures of this order may induce fatigue, leading to 
damage and an increased probability of breakdown 
after many rf pulses. 

4. Reverse bombardment model 

A fourth model [20] assumes that some ini- 
tial mechanism generates a local plasma (called a 
'plasma spot'), that by itself does not directly cause 



breakdown. Plasma spots, tiny sources of light, have 
been observed in DC and pulsed gaps without break- 
down. In this model, breakdown occurs when elec- 
trons emitted by the local plasma are returned to 
their source spots by the rf electric field. The en- 
ergy given to the source by these returning electrons 
is required to cause the plasma to grow and cause the 
actual breakdown. But as we will see later (Fig. 10), 
at least for 805 and 201 MHz in the absence of an 
axial magnetic field, electrons emitted at the high- 
est field location never come back to their source, 
no matter what their initial phase. This appears 
to be true for cavities in general, with the excep- 
tion of emission on the axis in pillbox cavities. Yet 
breakdown on the axis of pillbox cavities is rarely 
observed. 



5. Surface damage by heating 

As we have noted, breakdown gradients rise ap- 
proximately as the root of the rf frequencies; but this 
dependency does not continue at frequencies above 
10 GHz. Fatigue damage from cyclical surface heat- 
ing appears then to limit the gradients. This damage 
is worse at locations with maximum surface currents, 
where surface electric fields are usually low. But 
the damage, which causes cracks to form at grain 
boundaries can be so severe that breakdown is initi- 
ated in regions with both electric fields and surface 
current. Since this phenomena is only seen at such 
high frequencies, it is not expected to be a problem 
for neutrino factories or muon colliders. 



B. Observed Dependencies 

Although the above discussion suggest that ohmic 
heating initiates breakdown, this conclusion depends 
on a limited number of experiments; therefore one 
has to conclude that the initial mechanism that 
starts a breakdown is not fully understood yet. But 
there are some dependencies that appear fairly con- 
sistently: 

• Over a wide range of frequencies (0.2 to 
3 GHz), breakdown gradients are approxi- 
mately proportional to the square root of the rf 
frequency, and this dependency arises because 
of changes in the observed field enhancement 
factor ;^FN- 

• Breakdown appears dependent on the required 
energy to melt a given volume of the electrode 
material (see Fig. 7a). In the vacuum waveg- 
uide experiment case, this dependence is again 
caused by changes in /^fn, rather than in the 
local fields. 
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FIG. 7: (Color) Top: a) breakdown gradients vs the ap- 
proximate energy to melt a given volume of the material. 
Bottom: b) breakdown gradients vs the rf pulse length. 
(Data from Refs. 7, 18) 



• Breakdown occurs at lower gradients for long 
rf pulses than for short pulses (Fig. 7b), and 
there is some indication that this too arises 
from the pulse lengths influence on /3fn- 




FIG. 8: (Color) Breakdown gradients vs axial magnetic 
fields. Black line is dependency predicted by asperity 
twist model. Red lines are fit to the plotted Lab G [6] 
breakdown data. Blue lines are the calculated values for 
a 201 MHz cavity. Dotted lines are for Be surfaces. 



create bad asperities, leads directly to the observed 
dependencies with material properties and rf pulse 
length. 

In the case of the frequency dependency, since 
lower frequency cavities are usually larger, the ob- 
served dependency once again agrees with the ex- 
pectation. 



III. BREAKDOWN MODELS WITH 
EXTERNAL MAGNETIC FIELDS 

A. Published Breakdown Model with 
Magnetic Fields 



C. Conditioning 

In all cases discussed here, the cavities or waveg- 
uides were 'conditioned' prior to achieving the 
quoted gradients. After one or more breakdowns 
at one gradient, the cavity will subsequently with- 
stand a somewhat higher gradient. Typically, many 
hundred successive breakdowns are induced prior 
to the cavity reaching its 'final gradient'. When 
dark current distributions are studied as the gradi- 
ent is increased by this conditioning, then it is again 
found [21] that it is the enhancement /3fn that is 
decreasing and not that the local field is increasing. 

A reasonable assumption is that, whatever the ini- 
tial cause, the effect of a breakdown is to remove as- 
perities and consequently lowering /3fn- At the same 
time, however, the breakdowns, depending on the 
energy available, will create new asperities thus in- 
crasing /3fn- In a 'conditioned' cavity, these com- 
peting processes have reached an equilibrium. The 
rf energy available (stored or in a longer pulse), di- 
vided by the energy required to melt a given volume 
of metal, gives the total amount of melted metal. As- 
suming that more molten metal will be more likely to 



It has been proposed in the twist model [4] that 
the magnetic field dependence on breakdown arises 
from the torque forces on an asperity due to the in- 
flow of current feeding the field emission reacting to 
the external magnetic field; that is F (x I x B 
where E is the cavity electric field gradient and ap- 
proximately / oc £^^^ [4]. Thus for breakdown at a 
fixed force F, we expect 

£^breakdown OC B~^^^^ (8) 

In Fig. 8 the observed pillbox cavity breakdowns 
are plotted as a function of the external average 
magnetic field. The points plotted are those where 
both superconducting coils were powered so that the 
magnetic fields were relatively uniform over the cav- 
ity. The dependency predicted by this mechanism is 
shown by the solid black line; it is a poor fit at higher 
fields where the breakdown gradient falls much faster 
than predicted. 

B. Introduction to the Proposed Mechanism 

We propose a new model for breakdown with a 
magnetic field that is independent of the breakdown 
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FIG. 9: (Color) Proposed mechanism for breakdown 
with an external magnetic field. 



mechanism in the absence of magnetic fields (see 
Fig. 9). Breakdown occurs by this mechanism only 
if its breakdown gradient is lower than that from the 
case without a magnetic field. Its elements are: 

• 'Dark Current' electrons are field emitted from 
an asperity, accelerated by the rf fields, and 
impact another location in the cavity. In the 
absence of a magnetic field these impacts are 
spread over large areas and do no harm. 

• With sufficient magnetic field they are focused 
to small spots, where they can melt the surface 
producing local damage. If such damage is at 
a low gradient location there is no immediate 
breakdown, but the damage can accumulate 
until, for instance, a hole is made in a window. 

• If the electrons are focused onto a location 
with high surface rf gradient, then electrostatic 
forces will pull the molten metal out and away 
from the surface. This metal, as it leaves the 
now damaged location, will be exposed to field 
emitted electrons from the damaged area and 
will be vaporized and ionized, leading to a lo- 
cal plasma and subsequent breakdown. 

• For higher energy electrons, the melting will 
start deeper in the material where the ioniza- 
tion loss is greater, and expand to the surface. 
Thus, when the melting reaches the surface, 
significant quantities of molten metal can be 
sprayed onto other surfaces in the cavity [12]. 

Breakdown will be dependent on a) the Fowler- 
Nordheim field enhancement /3fn that determines 
the strength of the field emitted current, b) the local 
geometry of the asperity that will determine the ini- 
tial particle distributions and effects of space charge, 
and c) on the geometry and magnetic fields that fo- 
cus the electrons onto other locations. 



A program CAVEL [22] tracks particles from ar- 
bitrary positions on the walls of a cavity until they 
end on some other surface. The program uses SU- 
PERFISH [23] to determine the rf electric and mag- 
netic fields, and uses a map of external magnetic 
fields calculated for arbitrary coil dimensions and 
currents. Fig. 10 shows trajectories, for differing ini- 
tial rf phases, starting from the highest field location 
on an iris. Without an external magnetic field, none 
of the trajectories from the high field location come 
back to their common origin. Tracks emitted at a 
phase of 20° do hit the opposing iris at a high gra- 
dient location, but they are not focused there and 
are spread out over a significant distance. But with 
a sufficient external axial magnetic field, the tracks 
are either focused to the high gradient location on 
the opposite iris or returned to their source. 

There remains a small dependency of the arrival 
positions with phase, which arises from the com- 
bined effects of the perpendicular external and rf 
magnetic fields, but this dependency is small. If 
there were no other mechanism to spread out the 
electrons, then damage would appear as lines, which 
are not observed. In addition, damage would be 
much worse on the axis than at larger radii, also not 
observed. From this, one can assume that the dam- 
aging emission is dominated by field emission which 
is restricted to a limited range of phases (± ~ 20°), 
and that another mechanisms, such as space charge, 
spreads the beamlets out to a greater extent than 
this phase dependent effect. Figure 11 shows the en- 
ergies of the electrons on impact and indicates which 
phases are returned and which arrive on the next iris. 
Field emitted electrons (±25°) do not come back, 
but are focused to a local area on the opposite iris. 
Experimental observations [4] of the time structure 
of dark current show that this dark current is con- 
centrated at phases close to zero. It seems unlikely 
then that, even with an external axial magnetic field, 
significant currents of electrons will be returned to 
their source. 

The electron energies for the 805 MHz cavity with 
axial rf fields of 25 MV/m are approximately 1 MeV. 
For a 201 MHz cavity, and the same gradient, they 
are of the order of 4 MeV. At these energies, the 
electrons penetrate to significant depths in the Cu 
cavity walls. If Be is used, the penetration is even 
deeper. The relative surface heating and thus prob- 
ability of melting and damage depends on the frac- 
tion of energy deposited in a surface layer, taking 
into account thermal conduction away from the de- 
position. A full calculation would include the time 
dependence of both deposition and conduction. As 
an approximation, the temperature rise is estimated 
from the deposition in a depth corresponding to the 
thermal diffusion depth at the surface. 
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FIG. 10: (Color) Trajectories of electrons field emitted at different phases from the highest surface field location in 
an 805 MHz pillbox cavity with a) no external magnetic field, b) an axial field of 0.1 T, and c) an axial field of 1 T. 
The axial electric field is 25 MV/m. Phases are in degrees relative to the maximum. 
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FIG. 11: (Color) Energies of electrons on impact vs. 
their phase of emission. Red indicates electrons that 
returned, blue those that arrive on the next iris. Axial 
gradient is 17 MV/m. 



D. The Effects of Space Charge on the 
Transverse Distribution of Field Emitted 
Current 



Without an asperity and emission from a small 
area, the space charge forces give transverse mo- 
menta to emitted electrons causing the beamlet's ra- 
dius to increase. As the beamlet increases in radius 
and the electrons are accelerated, the space charge 
forces drop and it can be shown that the induced 
rms transverse momentum is (jp± ex vT. But if the 
electrons are emitted from the tip of an asperity then 
they will first be spread by the approximately spher- 
ically symmetric local electric fields, and the effect 
of the space charge is consequently modified. 

A simple simulation was performed (see Fig. 12). 
The initial electric fields were assumed to have 
strength PfnE, and exact spherical symmetry out 
to a distance X. Beyond this distance, the fields 
were assumed to be perpendicular to the average 
surface with strength E. No space charge effects are 
included in the spherically symmetric part. Radial 
space charge forces, inversely proportional to an av- 
erage radius, are assumed beyond the distance X. 
We found that with X = 0.2 //m, the transverse mo- 
menta could be approximated by ap± ex P , but with 
the power j = 0.3 instead of 0.5 as in the simple case 
without asperity. If X was smaller, then the best fit 
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FIG. 12: (Color) Schematic of approximate simulation of 
space charge effects on electrons emitted by an asperity. 




Gradient (MV/m) 



FIG. 13: (Color) The simulated final electron energy Se 
as a function of axial rf gradient for (red) a 805 MHz 
pillbox cavity, and (blue) a 201 MHz cavity. 



exponent j was found to increase in value. Since nei- 
ther the asperity height, nor its shape, are known, it 
is reasonable to treat the exponent j as an unknown 
that is fitted to the experimental data. In a better 
simulation, this fit would give us information on the 
asperity dimensions, but this model is too simple for 
this at this stage. At distances from the source large 
compared with X, space charge becomes negligible, 
but the transverse momentum is focused by the axial 
magnetic field, giving a beamlet with an rms radial 
size (Jr (X ^ where B is the axial magnetic field. 
The power per unit area W of the electron beamlet 
hitting the opposing surface is given by 



W 



(X 



(9) 



where £e is the final electron energy in (MeV) deter- 



mined from CAVEL [22] simulations, is plotted 
in Fig. 13 as a function of the axial rf gradient, for 1) 
805 MHz pillbox cavity [6], and 2) 201 MHz cavity 
designed for the MICE cooling experiment [24]. 



E. Fraction of Energy Deposited in the 
Thermal Diffusion Depth 

The thermal diffusion length S correspondin g to 
the rf field pulse duration r{s) is, S = 10~^VD r 
where D = is the thermal diffusion constant in 
(m), K is the thermal conductivity and Cs the spe- 
cific heat. It is assumed that all the energy deposited 
in this diffusion depth is spread uniformly within 
that depth S. The penetration depth d in (/im) of 
low energy electrons is approximately given by [25] 



d = 0.0267 



A 



1.67 



(10) 



where £e is the electron incident energy in (keV), p 
is the material density in (^^), and Z and A are 
the atomic number and atomic weight respectively. 

If £{x) is the energy of an electron with penetra- 
tion depth X and we define Q as the fraction of the 
electron energy deposited in the thermal diffusion 
length (5, then: 

• For electrons whose penetration depth is less 
than the diffusion depth {d < S): Q = 1. 

• For electrons whose penetration depth is 
greater than the diffusion depth, but not much 
greater {S < d < 10 x S) then: Q = 

Se-Se{d-5) 

• For electrons with penetration depth d > lOxS 



then: Q = 



F. Dependency of Local Temperature Rise 

The surface temperature rises as a fraction of the 
melting temperature is then 



AT 



W 



tQ 



S pCsT^ 



(11) 



where p is the density, Cs is the specific heat, Tm is 
the melting temperature, and r is the rf pulse length, 
taken to be 20 and 160 ps at 805 and 201 MHz 
respectively (r ex A^/^). Parameters used for Cu 
and Be are given in Table HI. 



G. Fit to Field Dependence Data from Lab-G 
Pillbox and Predictions for 201 MHz 

The above is a very approximate analysis. A full 
simulation of the problem is being pursued. The 



simulations were done with uniform magnetic field; 
tracks were only simulated from the single maximum 
gradient location; electron impacts were assumed at 
90° to the surface; the thermal diffusion calcula- 
tion ignored the rise time shape and used an ap- 
proximate calculation; both current scale and space 
charge strength were normalized to fit data. How- 
ever, if this simple approach can qualitatively fit the 
data, it should allow a qualitative extrapolation to 
201 MHz and other materials. A fuller simulation 
should provide more quantitative results. The red 
curved line in Fig. 8 shows the fit to the 805 MHz 
Cu cavity in Lab G breakdown data [6]; the fitted 
value of the current exponent was j = 0.35. More 
data [26] has been taken since the cavity had suf- 
fered significant damage. These data (not shown) 
lie at somewhat lower breakdown gradients and pre- 
sumably correspond to a now higher value of the 
Fowler-Nordheim /?fn- But since /3fn, has not been 
redetermined from dark current measurements of the 
damaged cavity, this analysis with that data is not 
possible. The black line in the plot shows the depen- 
dence predicted by the twist model [4], which does 
not fit the data well. The dashed lines show the cal- 
culated breakdown limits for Be. The horizontal red 
line indicates the gradient limit from the assumed 
model without magnetic fields, assuming the local 
field limit to be 7 GV/m. 



H. Extrapolation to 201 MHz 

Without an external magnetic field, breakdown 
gradients have been observed to follow approxi- 
mately a v7 behavior (see. Sec. II). Under the 
assumptions used here [17], the local field at break- 
down is independent of the frequency; this implies 
that /3fn decreases approximately as Using this 
assumption, and using the appropriately modified 
diffusion depth, the predicted breakdown limits at 
201 MHz, for Cu (solid blue line) and Be (dashed 
blue line) are shown in Fig. 8. The predicted 
201 MHz breakdown gradients are seen to be a fac- 
tor of 2 to 2.5 below those for 805 MHz. This factor 
comes primarily from the higher expected /3fn, but 
also from the longer rf pulse duration and thus longer 
time to heat and melt the surfaces. At the higher 
gradients specified for Neutrino Factory phase ro- 
tation and cooling (15 MV/m at B< 3 T) are well 
above this prediction. 



I. Other Experimental Results Consistent 
with this Analysis 

The pillbox cavity used in the above experimental 
study has, more recently, been tested [27] with one of 
its two Be window replaced by a flat Cu plate with an 



TABLE III: Atomic and Nuclear Properties of Cu and Be 
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Cs 
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^(805) 


^(201) 




J- m 








(W/cm-"C) 


(g/cm^) 


(J/g-°c) 


(cm^/s) 


(/xm) 


(/im) 




Cu 


29 


63.5 


4.01 


8.96 


0.385 


1.16 


48 


186 


1085 


Be 


4 


9.0 


2.18 


1.85 


1.825 


0.81 


40 


155 


1287 



easily replaceable central 'button'. The button had 
dimensions such that the local field on the tip was 1.7 
times that on the outer Cu 'iris'. The intent was to 
allow a study of magnetic field dependent breakdown 
as a function of materials, the assumption being that 
breakdown would occur at the high gradient on the 
tip of the button. The cavity was found to oper- 
ate with gradients on the button significantly higher 
(by approximately a factor of 1.7) than had been 
observed without the button. It was noted however 
that this breakdown occurred with gradients on the 
iris, and on the fiat support plate, that were essen- 
tially the same as those present without the button. 
This suggested that breakdown was not occurring on 
the button, but rather at other locations that did not 
have the 1.7 times field enhancement: on the irises 
and/or the fiat button support plate. When the cav- 
ity was later disassembled it was indeed found that 
was there little damage on the button, where the 
field was maximum. But there was significant dam- 
age in a distinct band from 3 to 6 cm on the TiN 
coated Cu support plate (see Fig. 14a). 

These observations are consistent with the break- 
down mechanism described in this paper. SUPER- 
FISH [23] calculations showed that the fields on the 
Be window opposite the button and support plate 
were maximal just in the band 3 to 6 cm, where the 
damage was observed (see Fig. 14b). In this model, 
emission from the Cu button, falling on the Be is 
less liable to cause damage and breakdown, whereas 
the emission from the Be focused onto the Cu plate 
should cause damage just where the gradients were 
maximal on the Be. 



IV. POSSIBLE SOLUTIONS, OR 
REDUCTIONS, OF THESE PROBLEMS 

There are several possible approaches to these 
breakdown problems in a Muon Collider or Neutrino 
Factory. 

• Redesign the phase rotation and cooling chan- 
nels to use lower rf fields. This approach would 
clearly hurt performance, and, in addition, 
risks a slow deterioration of performance as 
occasional breakdown events continue to spray 
Cu around the cavity - as observed in the pill- 
box tests. 

• Use cooling lattices with high pressure hydro- 
gen gas in the rf cavities. No degradation of 
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FIG. 14: (Color) Top: a) Button support plate showing 
damage band between 3 and 6 cm radius. Bottom: b) 
Surface fields on Be (red) and Cu (blue) vs radius, with 
bands showing where damage was concentrated. (Data 
from Ref. 27) 



rf performance has been observed in a small 
805 MHz test cavity with axial magnetic fields 
and rf gradients similar to those in vacuum 
cavities. If the loading from beam induced 
electrons is not a problem, or is slowed by 
introducing gas impurities, then this solution 
should offer no loss of performance in early 
cooling stages. But to cool to very low emit- 
tances, it will probably require lattices with 
lower Courant Snyder [5 at the absorber than 



-5 5 

Length (cm) 

FIG. 15: (Color) Open cavity with alternated solenoid 
coils in irises. 

can be achieved in the rf. In this case, the ad- 
dition of hydrogen gas at the higher ^s in the 
rf would cause unacceptable emittance growth. 
So for later cooling, high pressure hydrogen gas 
is probably not a solution. 

• Build cavities with exceptionally good surfaces 
so that the /3fn is sufficiently low initially that 
no breakdown occurs. With Atomic Layer De- 
position (ALD) this may be a realistic op- 
tion [28]. The fear would be that a single 
breakdown spoils the surface in such a way 
that there will be a cause of further breakdown 
and a conditioning that approaches the same 
gradient limits seen in a conventional cavity. 

• Design lattices with magnetic field shielded 
from the rf. The above prediction suggests 
that so long as the field is less than about 
0.2 T, no adverse effects will be observed. At- 
tempts to design such lattices have, unfortu- 
nately, shown significantly worse performance. 

• Design lattices using multi-cell open cavities, 
with alternating current coils in their irises 
(Fig. 15). In this case, as in the original 
multi-cell open cavity tests, the focused elec- 
trons would be directed to low field regions in 
the cavity and would thus not initiate break- 
downs. Nevertheless damage done to those lo- 
cations and molten Cu ejected from such dam- 
age could cause eventual deterioration of per- 
formance. In addition, the use of open, instead 
of pill-box cavities implies lower acceleration 
for a given surface field. 

• Damage might be eliminated if cavities were 
designed such that all high electric gradient 
surfaces were parallel to the magnetic fields 
(Fig. 16). This could provide 'magnetic insu- 
lation' [29]. Dark current electrons would be 
constrained to move within short distances of 
the surfaces, would gain little energy, would 
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FIG. 16: (Color) Top: a) Magnetic field lines from coils 
in a cavity lattice, together with cavity shape that fol- 
lows these field lines. Bottom: b) Energies of returning 
electrons as a function of their initial phase. 



cause no X-rays, and do no damage. Possible 
difficulties might be: a) cavities so designed 
will not give optimum acceleration for given 
surface fields, and b) multipactoring might oc- 
cur, now that the energies with which electrons 
do return to the surfaces are in the few hun- 
dred volt range where secondary emission is 
maximal. In addition, the use of open, instead 
of pill-box, cavities implies lower acceleration 
for given surface fields. 



V. EXPERIMENTS NEEDED TO STUDY 
THESE PROBLEMS 

Two critical experiments are already planned: 

1) The testing of the existing 201 MHz cavity in 
magnetic fields similar to those in current MICE 
and cooling designs. These experiments would use a 
MICE 'Coupling Coil' when it is available. 

2) Operating a high pressure hydrogen filled test 
cavity in a proton beam to study the possible break- 
down and rf losses due to the ionization of the gas 
by the beam. 

Whatever the results of these tests, further exper- 
iments will likely be needed to study the observed 
problems at 805 MHz and test possible solutions. If 
the tests of the 201 MHz cavity in magnetic fields 
show problems also at that frequency, then further 
805 MHz tests would be needed to explore solutions 
for these lower frequencies. Eventually, however. 
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FIG. 17: (Color) Simple pillbox cavity with mounting in 
two orientations within the lab G solenoid. 



tests of any proposed solutions would have to be 
done at 201 MHz. 

If with the proton beam the rf losses in the 
805 MHz test cavity are sufficiently low, then fur- 
ther tests of high pressure gas filled cavities will 
be needed. In particular, cavities must be tested 
with more stored energy and with a similar shape 
to those needed. Thin windows must be designed, 
safety problems must be addressed and beams with 
time structure and intensity nearer to those in the 
applications should be be employed. 



A. Experiments with a Simple Vacuum Pillbox 
Cavity 

The 'pillbox' cavities that have so far been tested 
have relatively complex shapes, and the simulations 
of their performance with fields in differing direc- 
tions are complex. It would thus be desirable to test 
a simple pillbox shaped cavity whose performance 
would be far easier to simulate. The cavity should 
be designed so that it can be mounted in the center 
of the of the Lab G magnet in at least two orien- 
tations: with its axis parallel to the solenoids and 
perpendicular to that angle. Ideally it should be 
possible to test it at other angles as well. A square 
cavity [30] (Fig. 17) would meet this requirements, 
would be easy to build and easy to arrange with 
coupling ports to the waveguide in either of two lo- 
cations to meet the two angle requirement. Such a 
simple cavity design would also be a good test vehi- 
cle for testing surface treatments including Atomic 
Layer Deposition (ALD). Since the design is so sim- 
ple, several versions could be made to compare their 
performances. 

As illustrated in Fig. 18, GAVEL simulations of 
a cavity with its axis perpendicular to the solenoid 
axis shows that the electrons are constrained to lie 
within a very small distance from the surface and 
gain little energy. Such a cavity would test whether 
multipactoring is a serious problem. In addition, if 
such a cavity were tested with the emission surface 
at a small angle to the field then the electrons could 
gain some energy, but would end up on the cylindri- 
cal outer surfaces where there would be no electric 
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FIG. 18: (Golor) A schematic illustration of an emitted 
electron's orbit when the magnetic field is parallel with 
the emitting surface. 




FIG. 19: (Golor) Single cell with magnetic insulation. 



field. This would test the configuration of the coil-in- 
iris solution without the shape modified to achieve 
magnetic insulation. Surface damage might be ex- 
pected, but it should not cause breakdown. 



A Single Cell Experiment with 'Magnetic 
Insulation' 



In this experiment (Fig. 19) superconducting coils 
would be mounted on either side of a cavity whose 
shape is such that the magneric fields are strictly 
parallel with the high gradient surfaces. The two 
coils close to the open pipe are powered in opposite 
directions. A second pair of coils is mounted outside 
the inner coils and powered with currents in the op- 
posite direction to those they surround. The coils 
would be separately powered so that the sensitivity 
to deviations from the insulated condition can be 
studied. 

The rf cavity would be operated at liquid nitro- 
gen temperatures in order to minimize the radia- 
tion falling on the superconducting coils that will be 



very close to its surfaces. There would be a vac- 
uum both outside the rf cavity (for thermal insula- 
tion) and inside the cavity, although the quality of 
the latter should be higher, and would be separately 
monitored. Super-insulation and an outer liquid ni- 
trogen shield are needed but not shown. Tuning of 
the cavity would be provided by axially squeezing 
or stretching the cavity. A variant of this experi- 
ment would have Be windows, making it into more 
of a pillbox design and raising the accelerating field. 
This would not provide such complete magnetic in- 
sulation but, because of the special properties of the 
Be might still have acceptable performance. 

C. A Multi-purpose Test Stand for these 
Experiments 

It is proposed that this and the following exper- 
iments would be carried out on a multi-experiment 
test stand. To allow ease of assembly, mounting 
of instrumentation, and making changes, all con- 
nections (rf, cryogens, magnet and instrumentation 
leads etc.) would be brought in through a single sup- 
port plate (on the right side in the illustration). The 
vacuum container would be in the form of a dome 
that joins to the support plate with a single flange, 
and can thus be easily removed without disturbing 
the connections. 



D. A Multi-cell Experiment with 'Magnetic 
Insulation' 




FIG. 20: (Color) Multi-cell cavity with magnetic insula- 
tion. 

The next experiment would be of a multi-cell mag- 
netically insulated cavity (Fig. 20). This would test 
the magnetic insulation in a geometry similar to that 
required in the 201 MHz acceleration for phase rota- 
tion or early cooling. In addition, it would address 
the problem of joining cavities inside the bore of the 
solenoids. This joint does not need to carry much rf 
current, nor does it need to hold high vacuum (since 
there is vacuum on both sides). The method requir- 
ing simple axial pressure is shown in the figure, but 
other ideas may be explored. 



It should be noted that this experiment, and the 
previous one, test a model of the 201 MHz lattic^| 
that would be used in the phase rotation. Neutrino 
Factory cooling, and the early cooling in a Muon 
Cohider. 



E. Testing of Components for later 6D Cooling 

Later 6D cooling will require compact higher field 
(10-15 T) solenoids in order to focus the beams to 
lower P's and thus cool to lower emittances. Test 
coils could conveniently be tested on the proposed 
test stand, since it would allow for easy mounting, 
cooling and testing with easily accessible instrumen- 
tation. Later, these coils would be used in conjunc- 
tion with 805 MHz cavities (see Fig. 21). The first 




FIG. 21: (Color) Lattice for later cooling to lower emit- 
tances. 

stage of the 10-20 T solenoid development could ap- 
propriately be tested in the same test stand dis- 
cussed above. Subsequently, the rf cavities would 
be included. The final step would be to add the re- 
quired hydrogen absorber. The details of the lattice, 
whether it is a snake [10] or Guggenheim [31] con- 
figuration for instance, would not change the nature 
of the components needing test. Figure 21 is thus 
not meant as a definitive design of the required lat- 
tice, but rather as a test of the kind of components 
needed; the cavity shapes shown are not those with 
true magnetic insulation. The true shapes with their 
coils are yet to be determined. Figure 21 is thus in- 
tended only as an illustration of the direction of the 
required R& D. 



F. Testing of Magnetic Field Solutions at 
201 MHz 

The experiments already planned with the 
201 MHz cavity in several Tesla fields could show 
satisfactory performance. If however this is not the 
case, then whatever solutions are proposed on the 



basis of the 805 MHz tests will need to be demon- 
strated at 201 MHz. 



VI. CONCLUSIONS 

The main results and predictions of the model pre- 
sented in this paper can be summarized as: 

• A review of models for breakdown without ex- 
ternal magnetic fields suggests that the initia- 
tion of breakdown is best explaned by ohmic 
heating due to field emission current at asper- 
ities. The dependence of breakdown on fre- 
quency, pulse length and cavity materials can 
be explained by the competing processes of as- 
perity destruction and the creation of new as- 
perities. 

• We have proposed a model for damage in rf 
cavities operated in significant axial magnetic 
fields. The model fits the existing data reason- 
ably well. 

• The model also fits some otherwise surprising 



results from a pillbox cavity with a Be window 

facing a Cu plate with a central button. 
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• The model predicts relatively low gradient 
breakdowns at 201 MHz in magnetic fields. 
These predicted breakdowns occur at signifi- 
cantly lower gradients than the operating gra- 
dients specified for the phase rotation and ini- 
tial cooling in the ISS [1] Neutrino Factory and 
in a Muon Collider. 

• Methods to address these problems are dis- 
cussed, including the use of 'magnetically in- 
sulated rf. An experimental program to study 
this concept is outlined. 
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